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We calculate the imaginary part of the induced pseudoscalar form factor of the nucleon 
Gp(t) in the framework of two-loop heavy baryon chiral perturbation theory. The effect of 
the calculated three-pion continuum on the pseudoscalar constant gp = (m IM /2M)Gp(t = 
— 0.877m 2 ) measurable in ordinary muon capture /j,~p — > v^n turns out to be negligibly 
small. Possible contributions from counterterms at twodoop order are numerically smaller 
than the uncertainty of the dominant pion-pole term proportional to the pion-nucleon 
coupling constant g n N = 13.2±0.2. We conclude that a sufficiently accurate representation 
of the induced pseudoscalar form factor of the nucleon at low momentum transfers t is 
given by the sum of the pion-pole term and the Adler-Dothan-Wolfenstein term: Gp{t) = 
Ag^NMU/iml - t) - 2g A M 2 (r 2 A }/3, with (r A ) = (0.44 ± 0.02) fm 2 the axial mean square 
radius of the nucleon. 
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The structure of the nucleon as probed by charged weak currents is encoded in two form 
factors, the axial and the pseudoscalar ones. To be specific consider the matrix element of the 
isovector axial current between nucleon states: 



{N(p + k)\qYl 5 T a q\N(p)) = u(p + k) 



l V G A (t) + ^G P (t) 
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where t = k 2 denotes the Lorentz-invariant squared momentum transfer and u(p) stands for 
a Dirac-spinor. M = 938.92 MeV is the (average) nucleon mass. The form in eq.(l) follows 
from Lorentz-covariance, isospin conservation and the discrete symmetries C, P and T. G^(f) 
is called the axial form factor of the nucleon and Gp(t) the induced pseudoscalar form factor 
of the nucleon. While experimentally much attention has been focussed on the first one, the 
latter is generally believed to be well understood in terms of pion-pole dominance as indicated 
from ordinary muon capture experiments fi~p — * v^n (see e.g. refs.[l, 2, 3]). The pseudoscalar 
coupling constant gp as measured in ordinary muon capture is denned via: 
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with tft the Lorentz-invariant squared momentum transfer if the proton and muon are initially 
at rest. m M = 105.66 MeV is the muon mass, M p = 938.27 MeV the proton mass, M n = 
939.57 MeV the neutron mass and = 139.57 MeV the charged pion mass. 

Chiral perturbation theory allows to calculate systematically the corrections to the dominant 
pion-pole term in Gp(t) (see eq.(7) below). At one-loop order this correction is uniquely 
expressed in terms of the mean square axial radius of the nucleon (r\) by making use of 
the chiral Ward identity of QCD [4]. Exactly the same term, derived originally by Adler, 
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Dothan and Wolfenstein [5] , is also found in the small scale expansion of ref . [6] where additional 
diagrams with intermediate A(1232)-isobars contribute. While the one-loop prediction for the 
pseudoscalar coupling constant gp = 8.4 ± 0.2 [4, 6] is consistent with the earlier result of the 
Saclay experiment gp = 8.7± 1.9 [2] a reanalysis [1] of that experiment using the modern world 
average of the muon mean life time gives the enhanced value gp = 10.6 ± 2.7. For further 
details on that and the conflicting situation concerning radiative muon capture fi~p — > u^ny 
see the recent review of Gorringe and Fearing [1] and also ref. [3]. 
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Fig.l: Two-loop diagrams contributing to the imaginary part of the induced pseudoscalar form 
factor of the nucleon Gp(t). Dashed and solid lines denote pions and nucleons, respectively. 
The wiggly line symbolizes the external isovector axial source. The combinatoric factor of the 
first four diagrams is 1/6 and the next four graphs have the combinatoric factor 1/2. The last 
two diagrams scale as g\ whereas the other eight graphs scale as g\. 

The purpose of the present short paper is to investigate the two-loop corrections to the 
induced pseudoscalar form factor Gp(t) in order to clarify whether these could affect (numer- 
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ically) the theoretical prediction for g P used so far. The essentially new feature at two- loop 
order is a nonvanishing imaginary part lmGp(t) for t > 9ml which originates from the (direct 
and indirect) coupling of the isovector axial current to the three-pion intermediate state. The 
pertinent ten topologically distinct two-loop diagrams generated by leading order vertices of the 
effective chiral Lagrangian + C^jj are shown in Fig. 1. The Feynman rules for the relevant 
interaction vertices can be found in appendix A of ref. [7]. 

Let us now turn to the evaluation of the imaginary part lmGp(t) from the two- loop diagrams 
shown in Fig. 1. Application of the Cutkosky cutting rules gives the spectral function lmGp(t) 
as an integral of the product of axialsource — > 3n and 3tt — > NN transition amplitude over the 
Lorentz-invariant three-pion phase space. Some details about these techniques can be found in 
refs.[8, 9] where the same method has been used to calculate the (two-loop) spectral functions 
of the isoscalar electromagnetic nucleon form factors and the 37r-exchange nucleon-nucleon 
potential. The choices e-k — and e v = k u for the polarization vector e u of the external isovector 
axial source allow a separation of ImGU(i) and Im[GU(i) + tGp(t)/4M 2 ]. Alternatively, one 
can use projection operators and tensorial integrals over the 37r-phase space can be reduced 
to scalar ones (see eq.(13) in ref. [8]). The pertinent three-body phase space integrals are most 
conveniently performed in the three-pion center-of-mass frame. The corresponding on-mass- 
shell four-momenta of the three pions read in this frame: k\ = (ui,ki), k 2 = (oo 2 ,k 2 ) an d 
fcg = (y/t — uji — ivi, —hi — k 2 ). The mass-shell condition k% = m\ determines the cosine of the 
angle between k\ and k 2 (called z) as: 
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zk x k 2 = uj x uj 2 - + w 2 ) + -(t + ml) , k 1>2 = s]u\ 2 - m\ . (3) 

The ten diagrams in Fig. 1 fall into two classes. The first eight diagrams carrying the common 
prefactor QaI ft give rise altogether to the following contribution to the imaginary part of the 
induced pseudoscalar form factor of the nucleon: 



t n (i), A _jhNM f , , f 1 ml Amj-ml mj(3ml - 1) 1m\u x u 2 zk 2 \ 

lmGp {t) ~WKr i 1 2 lT^"l2(t-m^ + ^o^ + (t-miy + t(t-ml)kj- 

Z 2 <1 

(4) 

Here, f n = 92.4 MeV denotes the pion decay constant and we have employed the Goldberger- 
Treiman relation: g n Nfn = 9aM. The inequality z 2 < 1 determines the kinematically allowed 
region in the Wi^-plane (which is bounded by a cubic curve) together with the obvious kine- 
matical constraints m n < u\ t2 < \fi — 1m v and Im^ < u\ + uj 2 < \fi — m n . Furthermore, 
one derives from the last two diagrams in Fig. 1 which are proportional to g\j f£ the following 
contribution to the imaginary part lmGp(t): 

g nN Mg\ r f 2 1— / k 2 \ aiccos(-z) 
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In the chiral limit = the total (two-loop) spectral function lmGp(t) shows a simple linear 
t-dependence of the form: 



t n u\\ 4M % n (A\ g * NMt 
ImG P (t)| m7r= o = — \m.G A {t)\ miT=Q 



t " v yim "- u 9(87r/ u ) 3 



2 / 047T^\ 



64vr 2 M St 

[ -W 2 - (6) 



3 



The first part of this equation follows from the fact that the combination + tG P (t)/AM 2 

is the form factor of the divergence of the isovector axial current, which vanishes in the chiral 
limit uin = (QCD chiral Ward identity). The (two-loop) result for lmGA(t)\m n =o has been 
taken over from eq.(27) in ref. [8]. 

In Fig. 2 we show by the full line the total imaginary part lmGp(t) calculated from eqs.(4,5) 
after division by a factor t. The horizontal dashed line in Fig. 2 indicates the asymptotic 
behavior of lmG P (t)/t for t — > oo. 
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Fig. 2: The spectral function ImGp(t) of the induced pseudoscalar form factor of the nucleon 
divided by t. The horizontal dashed line indicates its asymptotic form obtained by taking the 
chiral limit = 0. 

With the help of the spectral function lmGp(t) the complete two- loop representation of the 
induced pseudoscalar form factor of the nucleon can be written as: 
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The first two terms are the well-known pion-pole term and Adler-Dothan-Wolfenstein term 



[4, 5]. The parameters g-nNif-n 



TO, 



(r^) etc. are to be understood as the physical ones in- 



cluding their individual one- and two-loop chiral corrections. These (not explicitly calculated) 
two-loop renormalization effects come along with the real parts of the diagrammatic amplitudes. 
Note that the dispersion integral in eq.(7) requires two subtractions because of the asymptotic 
linear growth of the imaginary part lmGp(t). The third term in eq.(7) involving the two dimen- 
sionless low-energy constants Co and £i subsumes all polynomial contributions which arise from 
(tadpole-type) loop diagrams and possible chiral-invariant counterterms (beyond renormalizing 
the Adler-Dothan-Wolfenstein term). The prefactor of this term is chosen such that the negative 
mass dimension of the counterterm coupling strength is accounted for by appropriate powers of 
the chiral symmetry breaking scale A x = 2\f2-nf n ~ 0.82 GeV. Based on naturalness arguments 
one expects that the dimensionless low-energy constants Co,i are of order one. Indeed the same 



4 



considerations applied to the Adler-Dothan-Wolfenstein term give for the analogous dimension- 
less low-energy constant at one-loop order: — (2-Kf 7T ) 2 g A (r A )/3 = —gAi^Kf^jMA) 2 ~ —1.63. 
Here, we have inserted the value of the axial dipole mass Ma = (1.03 ± 0.02) GeV as extracted 
in ref.[10] from (quasi)-elastic neutrino and antineutrino scattering experiments. 

Let us now turn to numerical results. From the twice-subtracted dispersion integral in eq.(7) 
one gets a tiny contribution to the pseudoscalar coupling constant of 5g P ~ —1 ■ 10~ 5 . The 
extreme smallness of this number comes partly from the proximity of the chosen subtraction 
point to — to tp — — 0.502 m 2 . Nevertheless, when varying the subtraction point (via the 
substitution t 2 /t' 2 — > (t — t ) 2 /(t' — t ) 2 in eq.(7)) in the broad range —24ml < to < 9m 2 the 
contribution of the dispersion integral to the pseudoscalar coupling constant stays in magnitude 
smaller than one percent, > 5g P > — 1 ■ 10~ 2 . Note that a change of the subtraction point t 
is equivalent to changes of the low-energy constants Co,i parameterizing the polynomial piece 
in eq.(7). Since we cannot accurately determine the coefficients £o,i we turn here the argument 
around and ask only for some upper bound. For example, in order to cause a correction of 
Sgp = 0.09, corresponding to a relative 1% change of g^N in the pion-pole term, the relation 
2Co — Ci — 21 must hold. Such values of Co,i exceed the expectation from naturalness already 
by one order of magnitude. 

The delta-nucleon mass-splitting A = 293 MeV introduces another small scale to the prob- 
lem. The systematic power scheme inherent to the small scale expansion of refs.[3, 6] ensures 
however that the contribution to Gp{t) from two-loop diagrams with intermediate A(1232)- 
excitations is a homogeneous function of degree one in the three variables (t, m 2 , A 2 ). Conse- 
quently, possible negative powers of A get always overcompensated by two more powers of the 
numerically smaller scales m w and/or \f\t^\ in Gp{t^). 

One may therefore conclude that all two-loop corrections to the pseudoscalar coupling con- 
stant gp are numerically insignificant. A sufficiently accurate representation of Gp(t) at low mo- 
mentum transfers t is given by the sum of the pion-pole term and the Adler-Dothan-Wolfenstein 
term. Using for the 7rA^-coupling constant g nN = 13.2 ± 0.2 [11, 12] which is consistent with 
recent results from 7rA^-dispersion relation analyses [12] and (r A ) = 12/ M\ = (0.44 ± 0.02) fm 2 
[10] for the axial mean square radius one gets in this case: 

g P = 8.3 ±0.2. (8) 

The major theoretical uncertainty of gp comes obviously from the 7rA^-coupling constant g w N 
entering the dominant pion-pole term. 

Let us finally consider the form factor of the nucleon related to the divergence of the isovector 
axial current: 

The prefactor on the right hand side expresses the vanishing of this form factor in the chiral 
limit m n = as well as the presence of a pion-pole contribution. The imaginary part Imi5(t) 
(at two- loop order) can be easily constructed from the expressions for lmGp(t) given here in 
eqs.(4,5) as well as from the formula for lmG A (t) written in eq.(26) of ref.[8]. As a further 
non-trivial result we give here only the spectral function 1m D(t) in the chiral limit: 
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which may be useful for some quick order of magnitude estimates. 

In summary we have calculated in this work the imaginary part of the induced pseudoscalar 
form factor of the nucleon ImGp(t) at two-loop order in heavy baryon chiral perturbation 
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theory. Two-loop corrections to the pseudoscalar coupling constant gp measurable in ordinary 
muon capture fi~p — > u^n are numerically unimportant in view of the present uncertainty of 
the pion-nucleon coupling constant g n N- 

I thank T. Hemmert for useful discussions. 
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